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ABSTRACT 

We present a detailed spectral study (photoionization modelling and vari- 
ability) of the "Big Dipper" 4U 1624—490 based on a Chandra-High Energy 
Transmission Gratings Spectrometer (HETGS) observation over the ~ 76 ks bi- 
nary orbit of 4U 1624—490 . While the continuum spectrum can be modeled 
using a blackbody plus power-law, a slightly better fit is obtained using a single 
r = 2.25 power-law partially (71%) covered by a local absorber of column den- 
sity A^'h, Local = 8-l-a6 ^ lO^^cm"^. The data show a possible quasi-sinusoidal 
modulation with period 43^9^ ks that might be due to changes in local obscu- 
ration. Photoionization modeling with the xstar code and variability studies 
of the observed strong Fe XXV and Fe XXVI absorption lines point to a two- 
temperature plasma for their origin: a highly ionized component of ionization 
parameter ^hot ~ 10^'^ ergs cm s~^ (T ~ 3.0 x 10^ K) associated with an extended 
accretion disk corona of radius i? ~ 3 x 10^'^ cm, and a less ionized more variable 
component of ^ ~ 10^'^ ergs cm s~^ (T ~ 1.0 x 10^ K) and ^ ~ 10^'^ ergs cms" ^ 
(T ~ 0.9 X 10^ K) coincident with the accretion disk rim. We use this, with 
the observed Fe XXV and Fe XXVI absorption line variations (in wavelength, 
strength, and width) to construct a viewing geometry that is mapped to changes 
in plasma conditions over the 4U 1624—490 orbital period. 
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Subject headings: plasmas - accretion disks - X-rays: binaries, individual (4U 
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Introduction 



The "Big Dipper" 4U 1624—490, with its high energy (> 8keV) flaring between dips 
is one of the most 'extreme' of the dipping sources. Its nickna nie derives in part from the 
fact that it has one of the longest orbital periods (21 hr ft^ 76ks:lJones fc Watson 
dips that last ~ 2.7 hrs ^ llks (jWatson et al.lll985l : iBalucihska- Church et al. 



19891) and 
2o'ooir The 



prevailing idea for explaining this behavior is that the dips are likely due to obscuration 
by a raised accretion disk rim at the location where the infalling accretion stream impacts 
the disk. By accounting for the qu asi-steady X-ray dust scatt e ring halo, one ca n show that 
obscuration during dips is >90% (IXiang. Lee &: NowakI 120071 : llaria et al.l 120071 ). The dips 
exhibit broad, shallow shoulders, with a rapid ingress/egress to the deepest dip levels. A 
point-like soft X-ray source reasonably can be inferred from the rapid ingress/egress and dip 
variability. While the dipping seems to saturate at a minimum level, rapid (< 30 s), short 
duration, large amplitude spikes are seen. This behavior argues for the obscuring medium 
having a well-defined boundary, but a highly variable, erratic covering. 

19851) ■ vely with 



Since its discovery flWatson et al 



X-ray telescopes i ncluding EXOSAT ( Church fc Balucihska- Church 



2005 



19891 ). Eins tein ([Christian fc Swank 



BeppoSAX fjBalucihska- Church et al 



1997), ASCA flAngelini et al 



2000 



Smale. Church &: Balucihska- Church 



20061 ). and Chandra (llaria et al. 



2007 



11997 



19951). Ginaa (| Jones fc Watson 



Asai et a. 



, RX TE fiBalucihska-Church et al.ll2001 



Xiang. Lee fc NowakI l2007l : IWachter et al. 



2000), 



Lommen et al 



200lh. XMM -Newto n JParmar et al.ll2002l:lDfaz-Trigo et al 



20051). 



From these studies we know that 4U 1624—490 is highly absorbed, both extrinsically and in- 
trinsically. Models have been developed over the years wherein the spectrum consists of: large 
external absorption, a scattering halo, a point-like (i.e., rapidly abs orbed) blackbody, and an 
extended (i.e., partially absorbed) Accretion Disk Coro na (ADC) (jBalucihska- Church et al. 



2000 



2001; Smale. Church fc Balucihska- ChurchI 2001). Observations also have revealed a 



broad Fe Ka emission line (ISmale. Church fc Balucihska- ChurchI l200ll : iParmar et al.l |2002|) . 
It is therefore likely that we are viewing both the inner (point-like source, broad lines) and 
outer (dips, variable absorption) accretion disk regions in this highly inclined, high Edding- 
ton fraction, low-mass X-ray binary (LMXB). 

From the point of view of high resolution Chandra spectroscopy, 4U 1624—490 may be 
in many ways an alogous to the bright, highly inclined, variable, and absorbed microquasar 



GRS 1915+105 fiLee et all l2002l : iNeilsen fc Lee II2009I ). As such, there are key questions 
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concerning the physical environment of 4U 1624—490 that only high resolution X-ray spec- 
troscopy can address. Is the extended corona a photoionized a tmosphere in e quilib rium? 
Medium resolution XMM-Newton EPIC data provide some clues. iParmar et al.l (120021 ) iden- 
tified 'narrow' (<50eV) resonant absorption li nes from Fe xxvi Lyg, Fe xxvi He a, and a 
broad(er) 470 eV Fe Ka emission line (see also lSmale. Church fc Balucihska-Churchll200ll ). 
Like all X-ray sensitive CCD cameras, however, the XMM-Newton EPIC, with a resolution 
comparable to the ASCA SIS, is fundamentally limited in the extent to which it can resolv e 
such features. A recent analysis of these Chandra grating observations by llaria et al.l (120071 ) . 
while providing a better measurement of the ionized Fe lines given the better HEG resolution 
compared with EPIC in this spectral band, did not undertake the in-depth modeling neces- 
sary to describe the viewing geometry for the photoionized gas. Furthermore, based on the 
analysis presented in this paper, we offer an alternative scenario to the interpretation of laria 
et al. that the observed absorption line width can be attributed to bulk Comptonization. 

As a prelude to this paper, we studied the dust scattering halo properties of 4U 1624—490 
with this Chandra High Energy Transmission Grating Spectrometer (HETGS) observation 
in an attempt to first better understand the distribution of intervening line-of-sight material 
(size and spatial distribution of dust grains), which we foun d to be heterogeneously dis - 
tributed near the source and in the spiral arms of the Galaxy (jXiang. Lee fc Nowakll2007l ). 
Additionally, by measuring scattering delay times as a function of angular distance in com- 
parison to the 4U 1624—490 source light curve, we were able to obtain a much more accurate 
distance measurement of 15^2 e kpc to 4U 1624—490 . In this paper, we use the same Chan- 
(ira-HETGS observation to focus on the binary propertie s with an in-dep th analysis of the 
HETGS spectra. (These are the same data published by llaria et al.ll2007l .) Specifically, we 
explore the persistent-phase spectrum over the 4 U 1624—490 o rbital period in the context 
of §2.11 photoionization modeling using xstar (IKallmanl |2005| . Version 2.1kn7), and §2.21 
detailed variability studies of the detected ionized Fe XXV and Fe XXVI absorption lines and 
continuum. In particular, we study the orbital phase evolution of the detected absorption 
lines as a function of the time periods immediately before and after ( §2.2.ip . as well as near 
to and far from the dipping events ( §2.2.2p . in addition to an analysis of the short-duration 
(^ 10 ks) absorption changes over the 4U 1624—490 ~ 21 hr orbital period. The results are 
summarized and discussed in §3] in the context of a simple two temperature absorber model, 
and we specifically focus on the quasi-sinusoidal modulation of these plasma components, as 
well as modulation of the continuum components, over the binary orbital period. 
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2. Observations and Data Analysis 

We observed 4 U 1624—490 with the Chandra High Energy Transmission Grating Spec- 



trometer (HETGS, ICanizares et al.ll2005l ) in timed graded mode beginning at 06:26:06 UT 



on 2004 June 04 (MJD 53160.26813, ObsID: 4 559), covering one binary orbit of ^76 ks . 



Given that 4U 1624-4 90 is bright (?^50 mCrab. ISmale. Church fc Balucihska-ChurchI 12001 



Diaz-Trigo et al.ll2006l ). we mitigate pileup and telemetry problems by using a one-half sub- 
array of 512 columns per CCD in order to reduce the frame-time to 1.7 s from the nominal 
3.2 s. Despite these efforts, the 0*^ order ACIS-S data are heavily piled-up, rendering diffi- 
cult the drawing of useful conclusions from the 0th order spectrum. Accordingly, this paper 
focuses only on spectral and light curve analysis of the HETGS grating spectra. 

We use CIAO 3.4 with caldb 3.3 for our data reduction efforts. In order to ensure 
the most accurate wavelength measurements possible for this observation, we calculate the 
exact zero order position using a code that fits the intersection of the dispersed MEG data 
with the zeroth-order "readout trace" on the ACIS-S3 chipQ. Using this technique it is 
expected the accuracy of the O*'' order position is better than 0.1 detector pixels, translating 
to a wavelength accuracy of 0.001 A (corresponding to 164 km s~^ at 6.8 keV) for both the 
High Energy Grating (HEG) and Medium Energy Grating (MEG) first orders. Positive and 
negative first order spectra for HEG and MEG were combined for analysis. 

For our analysis, we generate f^^-order time-averaged HEG and MEG spectra spanning 
the non-dip phases of the full 76 ks orbital period (resulting in ^Q7 ks of usable data) 
to study the properties of the persistent phase continuum spectrum and to look for faint 
lines ( §2.11) . We also generate phase-dependent spectra (labeled "a-g", and combinations 
thereof, as indicated by the light curve in Fig. [T]) to study the time evolution of the plasma 
via observed changes in the line profiles and fluxes ( §2.21) . Figure [1] also shows that our 
observation encompasses ~9.7 ks (~2.7 hrs; sum of dip phases 1-1-2-1-3) worth of occulted 
periods. The longest duration of an individual dip is seen to be ~3.5 ks, and the largest 
occultation is ~90%. Possibly because of the low energy absorption, the spectrum is seen 
to harden during dipping periods. We are unable to extract useful information from the 
dip-phase HETGS spectrum due to the low counts. 



^ ht tp : / /space . mit . edu / cxc / analysis /findzo / 
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2.1. The Average Persistent Phase Spectrum 



As a prelude to spectral analysis, we first determined pileup fraction (defined to be the 
fraction of frames that have detected events containing two or more events per frame per 
pixel) in the first order grating spectra. Since pileup is dependent on the count rate in each 
CCD pixel, we use the count rate in the HEG±1 and MEG±1 spectra for our estimates. We 
find that the maximal pileup fraction of HEG +1 and —1 order spectra is ~1.0% at ^2.5 A 
(^5 keV). Similarly MEG +1 and —1 order spectra are piled up at less than f5il.0% except 
for the range 2.5-3.5 A (3.5-5.0 keV) where the maximal pileup fraction reaches ~2.0%. 
Accordingly, we ignore regions of the spectrum with pileup >1%. As such, we fit only to 
1.2-7.3 A (1.7-10.0 keV) HEG spectra, and MEG spectra over the 1.5-2.5 A (1.7-3.5 keV) 
and 3.5-7.5 A (5.0-8.0 keV) range. The spectra are rebinned to satisfy a signal-to-noise 
ratio S/N > 5 and 0.005 A per bin for the HEG 1^* order 5/ A > 5 and 0.01 A per bi n for 
the MEG order, respectively. ISIS@ faouck fc Denicolalboool : iNoble fc Nowaklboosh was 
used as the fitting engine for our analysis. 



2.1.1. Broadband Continuum 



The 4U 1624—490 broadband continuum as determined from the HEG spectrum has 
already been discussed at length in our previous pap er focusing on the scattering halo as- 
sociated with this source (jXiang. Lee &: NowakI 120071 ). The inclusion of the MEG spectrum 
to our fitting here gives consistent results for a blackbody + power-law continuum modi- 
fied by absorption from neutral material (Ah). These parameters are listed in column 3 
of Table [H and a spectrum is shown in Fig. [2l In our previous work, however, we were 
concerne d by the fac t that our best-fit parameters were noticeably different from that pub- 
lished by llaria et al.l (120071 ) for the same data set. In particular, our best-fit photon index 



r = 1.14_Q'4g is slightly harder than the F = 2.0_q'55 noted by laria et al. A cursory check 
showed that the difference had very little to do with absorption models, e.g., F f» 1.4 for 
PHABS (IBalucihska-Church fc McCammon Ill992h. versus F 1.2 f or tb abs, which use the 
updated cross-sections and abundances of IWilms. Allen fc McCrayl fl2000l ). fitting tools (e.g., 
ISIS versus XSPEC), binning, or choice of fitted spectral regions. For completeness, we also 
looked into calibration effects and find that at least some of the differences had to do with 
a difference between caldb 3.2 (used by laria et al.) versus caldb 3.3 (this paper). A 
reduction of the 4U 1624—490 spectra using caldb 3.2 resulted in a slightly softer power- 
law (F 1.8), compared to our caldb 3.3 value of F 1.4 (phabs). Even though the 



^http: / /space, niit.edu/cxc/isis/ 
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Light Curve of 4U 1624-490 
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Fig. 1.— (A) First-order 1-10 keV HEG+MEG light curve of 4U 1624-490, and (B,C) 
associated hardness ratios that reveal spectral hardening during the dip periods. Throughout 
the paper, time zero corresponds to MJD 53160.3. 
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Wavelength (A) 

Fig. 2. — The best-fit persistent phase continuum and broad Fe emission hne model (red) 
of Table [1] modified by two photoionized absorption lines (see §2.1.31) overplotted on the 
time-averaged HEG 1*^* order spectrum of 4U 1624—490 over its ~ 76 ks orbital period. The 
blue line is the best-fit model excluding the Fe emission line and the ionized absorption lines. 
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CALDB 3.2 measure of photon index is closer to the T = 2.0!l^o'55 value quoted by laria et al. 
for this data set, a full agreement with the values of laria et al. is not reached even when 
trying to reproduce every step of their analysis process. We can only speculate that a com- 
bination of effects, including new reprocessing that accounts for afterglow events, and the 
careful treatment of pileup incorporated here may be responsible for the rest of the discrep- 
ancy. We also note that our ph oton index is slightly har d er th an the values of F = 2.0lo 8 
for a BeppoSAXohseivation bv iBalucihska- Church et al.l (120001 ) and F = 2.02 ± 0.12 for a 
XMM-A'^ewton observation by lParmar et al.l (120021 ). The difference is possibly caused by the 
different luminosity states (see the next paragraph) or by the limited Chandra-HETGS band 
coverage (e.g. 0.5-10 keV) and the severe absorption at the low energies coupled with limited 
counts at the high energies. 



Based upon our best-fit continuum, the 0.6-10 keV flux of 4U 1624—490 is ~ 1.07 x 
50 mCrab), corresponding to an unabsorbed 0.6-10 keV luminosity 



10 ^ erg cm ^ s ^ 



4.9 X 10^^ er gs s ^ (assuming a distanc e of 15.0 kpc as determined from our scatter- 



of 

ing halo analysis; IXiang. Lee &: Nowakl l2007r) . The luminosity is similar to that observed 
with XMM-Newton (jParmar et al.ll20021) bu t significantly low er than that using the RXTE 



(iSmale. Church fc Balucihska- Church 
I995I ) observations. As discussed by 



200lh and EXOSAT (IChurch fc Balucihska- Church 



Smale. Church fc Balucihska- ChurchI (1200 ll ). the lu- 



minosity during different observations can change significantly and thus the emission param- 
eters should not necessarily be expected to completely consistent. 

As an alternative to the continuum model of a blackbody plus power-law, we also 
considered a continuum consisting solely of a partially covered power-law. In this model, the 
soft excess that had been described by the blackbody is now modeled by the soft end of the 
uncovered power-law. Specifically, the model is given by 



X 



xe 



H, LocalO"ph(-E) j 



where A^h, los is the equivalent hydrogen column along the line-of- sight, A'h, Local is the 
equivalent hydrogen column local to the system, A is the normalization of the power-law 
continuum, and x is the partial covering fraction. We find that this model {x^/d.o.f = 
1413/1338) describes the persistent phase spectrum slightly better than the blackbody plus 
power-law model {x^/d.o.f = 1435/1338), and yields values of x = (71±3)% , F = 2.25±0.06 
and A'h, Local = 8.llo;g x 10^^ cm~^. Again, as discussed above, the limited Chandra-BETGS 
band coverage (e.g. 0.5-10 keV) and the limited counts at high energies do not allow us to 
distinguish between these two models for the time averaged spectra. Thus for some of the 
results that follow, e.g., the modeling of the broad emission line ( §2.1.2p and narrow ionized 
absorption lines ( §2.1.31) . the results are insensitive to which continuum model that we choose. 
As we elaborate upon in the discussion of short time scale changes ( §2.2.3p . however, the 
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partial covering model may more naturally explain spectral variations associated with the 
source's light curve. 



2.1.2. Broad Emission Line 



A significant and broad excess around 6.6 keV was found when we checked the fit 
residuals in both the HEG and MEG spectra. We also saw deficits around 6.7 keV and 
7.0 keV corresponding to narrow absorption lines. These are discussed in more detail in 
§2.1.31 and §2.21 The broad excess has be en identified with Fe K emission typically de- 



tected in LMXBs (1*^* significant detection: ISuzuki et al 



White et al. 1986 and Hirano et al. 



1987; 



Asai et al 



2000 



1984 with Tenma: with EXOSAT, 



for an ASCA study of 20 LMXBs, 
including 4U 1624— 490 ). I Asai et al.l (120001 ) postulate that the line is likely produced through 



the radiative recombination in a photoionized plasma and that its broad width is attributable 
to a confiuence of processes in t he ADC, including C ompton scattering and Doppler shifts 
arising from Keplerian motions. iParmar et al.l (120021 ) also studied the broad emission line 
from 4U 1624—490 with XMM-Newton , and found that the best-fit energy for the emission 
line depends on whether or not the Fe abs orption lines were also modeled. When the ab- 



sorption lines were included, iParmar et al.l found a best-fit emission energy of 6.58lo oJ keV 



in contrast to the 6.39lo;o4 keV value found when the absorption features were excluded. 
The HEG clearly shows strong Fe XXV and Fe XXVI absorption lines in our HEG spec- 
trum, which are therefore included in our modeling of the the broad emission feature. The 
continuum model is the blackbody plus power-law model of §2.1.1[ 

Like the persistent phase spectrum, the far-dip {a+e+f+g; ~39 ks) and near-dip {b+c+d; 
~29 ks) spectra of Fig. [3] reveal a broad emission line with a best-fit energy that is most 
consistent with He-like Fe XXV at ~ 6.7 keV (Table [1]). It is interesting to note that the 
peak energy for the broad emission feature in the near-dip phase spectrum points to even 
higher ionization H-like Fe XXVI and to a broader velocity width, although the uncertainties 
of the modeling do not allow us to make any strong claims to this effect. If real, it would be 
of interest to consider whe ther Compton sca ttering can be a mechanism for line broadening 
in LMXBs as suggested by lAsai et al.l (120001 ). Note that the highest values for the fitted line 
width, 0.5 keV, correspond to velocities that if associated with the accretion disk arise 
at radii > 100 GM/c^. Thus, the broad emission line is not necessarily cospatial with the 
blackbody component, if indeed a blackbody-|-powerlaw is the better description of the true, 
underlying physical description of the observed spectra. The presence of a broad line there- 
fore does not provide us with a means of distinguishing between the blackbody+powerlaw 
and partially covered powerlaw models for the continuum spectra, and could be consistent 
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with either modeL 



2.1.3. Photoionization modeling of the narrow ionized absorption lines 



The three most obvious absorption hnes that show up in the Chandra-HETGS spectra 
(Fe XXVI at ^ 7.0 keV/1.778 A Fe XXV at keV/1.850 A and Ca XX at ^4.1 keV/3.020 A; 
see Fig. |2]) have aheady been reported by laria et al, albeit based upon Gaussian fits to these 
data, as opposed to the modehng with photoionization codes as we present in this work. lon- 
i zed absorption hnes were initially discovered for 4U 1624—490 in XMM-NewtonEPlC data 
(jParmar et al.ll2002l ). Here, for the time-averaged spectra we go beyond previous studies 
by modeling the absorptio n spectrum of 4U 1624—490 using an analytic version of xs- 
TAR v2.11n2 (jKallmarul2005l ). For these xstar fits, the column density, ionization, turbulent 
velocity width, and velocity shift were allowed to vary. The abundance of elements from 
C to Fe were fixed at solar values (iGrevesse et al.lll996l ). although subsequent fits that al- 
lowed these parameters to vary gave abundance results similar to solar values within 90% 
confidence errors. If we adopt the interpretation of laria et al. and consider a one-zone 
plasma, the time-averaged spectrum over the course of the 4U 1624—490 orbital period 
can be modeled reasonably well with a hot plasma of log C, = 3.6 ± 0.2 ionizing a col- 

Be- 



656l^^°kms-i. 



umn of log iVn = 22.53to;o8 cm~^, with turbulent velocities of Vturb 
cause our variability analysis of §2.21 strongly points to a two-ionization zone model for the 
4U 1624—490 absorbers, we describe our model for the time-averaged spectrum accordingly. 
For this model, a combination of ionized absorber-1 (log^ = 4.3 ± 0.4; log A^h = 23.3 ± 0.2; 
vturh = 280l^^°kms-\ and v.^it = -607l^^^km s"^), and absorber-2 (log^ = 3.3 ± 0.2; 
logA^'n = 22.1 toicm"^; fturb < 174kms~^ and fshift = 213lj;58kms~"'^), provides a better 
fit. Interestingly, for either scenario the Ar XVIII component cannot be modeled well unless 
the abundance of Ar is set to f^^twice solar. Table [2] lists ions detected in our spectrum, 
and also the strong lines that are predicted by our xstar models to have equivalent widths 
Priori > 0.1 mA, including three strong lines of Fe XXVI, Fe XXV and Ca XX. 



2.2. Orbital Phase Evolution of the H-like and He-like Fe Absorption 

Given that the light curve (Fig. [1]) clearly indicates a blockage of the primary by the 
accretion stream and/or secondary, we next investigate relevant effects this might have on the 
observed ionized absorption. For this, we concentrate only on the evolution of the strongest 
absorption lines, i.e., Fe XXV and Fe XXVI, over the course of the 4U 1624—490 ~ 76 ks 
orbital period, as the source moves in and out of occultation. For this analysis, we use only 
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Table 1. Best-fit continuum and broad Fe emission line parameters 



Component 


Parameter 


persistent 


pre-dip 


post- dip 


near-dip 


far-dip 




x'/d.o.f 


1435/1338 


1003/992 


1205/1162 


1026/1066 


1248/1194 




Lx (erg s-i) 


4.9 X lO^'^ 


4.8 X 10^'' 


5.0 X 10^'' 


4.6 X 10^^ 


4.9 X 10^'' 


TBABS" 


Nh (10^2 cm~2) 


'■"-0.4, 


7 4+0.6 
'■^-0.4 


7 0+0.5 
'-S-0.5 


7 c+0.6 
'•^-0.5 


' -"-0.4 


Blackbody 


kT (keV) 


1 oQ+0.05 
J^--3»_0.04 


1 qo+0.08 


1 4Q+0.O6 

l-*'^-0.05 


1 40+O.O8 
l-^'^-0.06 


, qc: + 0.06 
^^•■^■^-0.04 




Norm (10^2) 


1 42+015 
1-^^-0.11 


, c-r,+0A6 
^•■^^-0.29 


, c:o+0.14 
l-'J'-'-0.12 


1 45+0.27 

, ,0+0.54 
1-1'^-1.17 


1 41+0.24 


Powerlaw 


r 


1 14+0.33 


0.651?:?^ 


-1 4-3+0.43 
l-*'^-0.55 






Norm (10^2 ph cm-^ kcV"^) 




r> rrH-8.3 

^- ' -2.6 


q 1+9.9 


= 1+8.8 
0.i_4.7 




Fe emission 


Efc (koV) 


6.76t0;?? 


6 57+0-19 


O.DO_0.08 


6 Q7+0-47 
"■^'-0.31 


6.68iH^7 


feature 




(1.83 A) 


(1.89A) 


(1.86A) 


(1.78A) 


(1.86A) 




Flux /pc (10-3 ph cm-2 s-l) 


^■='-0.3 


0.8«i 




1 9+0.6 
^■^-0.6 


0.8«-J 




Width a (keV) 


0.341°;}! 


0.24±«-9 


0.13tO;g? 


50+0-32 


18+0-07 
^•1^-0.05 




(10^ km s-i) 


(151^) 




(6tJ*) 


(22^1-1) 


(8t^) 



Note. — "tbabs is the absorption model described in ^2.1.11 Errors arc quoted at 90% confidence. 



Table 2. Ions predicted from xstar fit to time averaged spectrum 



Ion" 




Transition* 






obs 


rf 




model'' 


comments 


Mg XII 


Ly /3: 




» 3p(2p°) 


0.079 


7.106 


7.111 


6.75 


0.51 


w2 


weak 


Si XIV 


Ly a: 




<■ 2p(2p°) 


0.414 


6.182 


6.170 


5.45 


1.10 


wl 


weak 












6.186 


81.70 


4.10 


w2 


strong 


Si XIV 


Ly 13: 




» 3p(2p°) 


0.079 


5.217 


5.221 


13.70 


0.70 


w2 


weak 


S XVI 


Ly a: 




» 2p(2p°) 


0.413 


4.729 


4.720 


6.31 


0.98 


wl 


weak 












4.732 


83.00 


3.20 


w2 


strong 


S XVI 


Ly /3: 




> 3p(2p°) 


0.079 


3.991 


3.994 


14.20 


0.58 


w2 


weak 


Ar XVIII 


Ly a: 


ises,/2) - 


» 2p(2p°) 


0.412 


3.733 


3.736 


22.50 


0.79 


w2 


weak 


Ca XX 


Ly a: 


ises,/2) - 


» 2p(2p°) 


0.411 


3.020 


3.022 


12.00 


0.62 


w2 


weak 


Fe XXV 


He a: 




ls2p(ipf ) 


0.775 


1.850 


1.851 


113.00 


2.30 


w2 


strong 


Fe XXV 


He /3: 




ls3p(ip°) 


0.152 


1.573 


1.574 


31.00 


0.74 


w2 


weak 


Fe XXVI 


Ly a: 




* 2p(2p°) 


0.408 


1.780 


1.776 


103.00 


4.00 


wl 


strong 












1.781 


145.00 


2.10 


w2 


strong 


Fe XXVI 


Ly /3: 




* 3p(2p°) 


0.078 


1.503 


1.500 


12.70 


0.72 


wl 


weak 












1.504 


29.60 


0.65 


w2 


weak 


Fe XXVI 


Ly 7: 




> 4p(2p°) 


0.029 


1.438 


1.435 


5.14 


0.26 


wl 


weak 












1.439 


10.90 


0.26 


w2 


weak 


Fe XXVI 


Ly S: 




■ 5p(2p°) 


0.014 


1.404 


1.401 


2.65 


0.12 


wl 


weak 












1.405 


3.77 


0.13 


w2 


weak 



Note. — '^Oscillator Strength, ''Laboratary Wavelength, '^Observed wavelength, ^Optical Depth, f Equivalent Width, ''"wl" 
and "w2" refer to the absorber-1 and absorber-2 respectively. 
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the HEG given its higher resolution and throughput in this energy band. 

2.2.1. Pre-dip versus Post-dip phases Fe line changes 

A simple overlay of the pre-dip {a+b; ~20 ks) and post-dip {d+e+f+g; ^38 ks) phases 
reveal nearly identical spectra. Our analysis here uses the blackbody plus power-law contin- 
uum ( §2.ip and broad Fe emission line ( §2.1.2p to model the more global features. While the 
underlying components are similar to that used for the time averaged spectrum, the best-fit 
parameters differ slightly. To quantify this effect, we assess Fe line changes by adding Gaus- 
sians to model the Fe XXV and Fe XXVI absorption. The best-fit parameters (Table [3]) do 
not indicate significant changes in absorption between the pre-dip and post-dip phases, so 
for these fits we do not proceed with the more involved xstar modeling. 

2.2.2. Near-dip versus far-dip phases Fe line changes 

Next, we considered a comparison between the near-dip {b+c+d; ~ 29 ks) and far-dip: 
(a+e+f+g; ~ 39 ks) spectra following (initially) the procedure of §2.2. 1[ Based on simple 
Gaussian fits, distinct Fe line changes are observed in the line parameters, most noticeably 
in the Fe XXV velocity width, equivalent width, and line flux (Fig. [3] and Table |3]). Fig. [3] 
and Table [3] reveal distinct Fe absorption line changes with the Fe XXV:Fe XXVI flux ratios 
indicating noticeable plasma changes (temperature & ionization) between the near-dip and 
far- dip phases. 

In an attempt towards a more physically realistic picture explaining these changes, we 
proceed with more rigorous photoionization modeling of these spectra. For this, we fixed the 
continuum parameters specific to each phase, and varied the xstar parameters as before to 
model the Fe XXV and Fe XXVI absorption features. For these spectra, we find that a single 
ionization absorber gives an overall reasonable fit in a sense, but does not fully describe 
the Fe xxv absorption line. The fact that the relative Fe XXV and Fe XXVI velocity shifts 
are different in the far-dip and near-dip spectra (i.e.. Fig. [3]) suggests that a more realistic 
physical model should consider two ionized absorbers arising from different physical regions. 
A fit to these near-dip and far-dip spectra using such a two-component ionized plasma: one 
hot (~ 3.0 X 10^ K corresponding to log.^ ~ 4.2) and another slightly cooler (Ri 1.0 x 10^ K; 
log^ ~ 3.0), gives a better description (Fig. [3]). By adding this second component, the 
of the fit to the near-dip phases is reduced from 1047 for 1063 d.o.f (one-component ionized 
plasma) to 1033 for 1059 d.o.f (two-component ionized plasma). This yields an F-statistic 
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value of 3.6, which for these four additional degre es of freedom formally corresponds to an 
improved fit at > 99% confidence. (See, however, iProtassov et al.l |2002| for a discussion of 
the applicability of the F-test in situations such as those presented here.) We therefore view 
the two-component ionized plasma model as a more likely scenario for the 4U 1624—490 
X-ray emitting regions. 



2.2.3. Short duration changes 

Since there is clear evidence for absorber changes when comparing the far- and near-dip 
periods, we conclude our analysis of the spectral variability by dividing the total persistent 
spectra into 7 approximately equal duration parts of ~10 ks each (Figs. [1], [Hand [5]). Of 
interest is whether the blackbody-|-power-law or partial coverer-|-power-law components vary 
with orbital phase. We first test for such variations by fitting our blackbody plus power- 
law continuum model, including the broad emission line and narrow Fe absorption lines. 
Fig. [3]-(II, III) and Table O show that the unabsorbed fluxes of both the blackbody and 
power-law components, and to some extent the absorbed flux of the blackbody component, 
remain relatively steady over the 4U 1624—490 orbital period. The absorbed power-law 
component, however, is seen to vary more dramatically, and in a quasi-sinusoidal fashion 
(~95% confidence). 

We note that the quasi-sinusoidal variation of the absorbed blackbody component is 
less than that of the absorbed power-law, and perhaps exhibits an anti-correlation with the 
power-law. We therefore assess the statistical correlation between these parameters by gener- 
ating confidence contours for phase (g) (see Fig H] and E]) . These contours indicate a marked 
degree of dependence of the absorbed power-law and blackbody fluxes upon each other; 
therefore, in deriving a period for the modulation of the flux, we tie these two components 
together and fit a sine function to the points of Fig. [5]-(Ill). Specifically, the power-law and 
blackbody fluxes are tied together, but the power-law flux is allowed a 1/2 period (delayed) 
phase shift. At 90% confidence, we derive a period T = 43^1^ ks (~ 12 hrs), slightly more 
than half the orbital period. 

Next, we explore the time dependence of the partial covering model (see the description 
in §2.1.11) by fitting simultaneously the seven phased spectra. We wish to determine whether 
the change of the absorbed power-law flux is produced by variation of the local cold gas. 
Here we tie together all fit parameters (including power-law normalization and slope), but 
allow the partial covering fraction and local absorption, i.e, the cold absorption along the 
line-of-sight, to vary with orbital phase. We find that this model describes the seven phased 
spectra well, giving of 5971 for 5896 degrees of freedom and yields values of iVn = 
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Far Dip 




01.6 1.7 1.8 o 1.9 

Wavelength (A) 

Fig. 3. — Fluxed spectra (black) and best-fit two-component ionization model (red) of §2.2.21 
during the far-dip (top) and near-dip (bottom) 4U 1624—490 phases. Represented in light 
blue is a single-component ionization model that clearly fits the data less well, in dark blue is 
the model excluding the ionized absorption lines. Vertical green lines indicate the laboratory 
wavelengths of Fe xxv (1.850A) and Fe xxvi (1.778A). 
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Table 3. Best-fit Fe xxviand Fe xxv absorption lines az measured in the different phases. 



Orbital 


Wavelength 


Width 


EW 


Line 


Wavelength 


Width 


EW 


Line 


phase 


(A)(keV) 


(mA)(km/s) 


(mA) 


Flux 


(A)(keV) 


(mA)(km/s) 


(mA) 


Flux 






Fe XXVI 








Fe XXV 






nondip 
predip 
postdip 
neardip 
fardip 


1.778tH°02(6-97) 
1-7801°:°m(6-97) 
1.7751H;3'2(6-99) 
1.778tHo?(6-97) 
l-778iHo3(6-97) 


5i^(770lfJ) 

< 7.0{< 1176) 
4ti(758t«|) 

< 3.8(< 635) 
6«(1094t«?) 


7 9+1.3 
'■^-1.5 
c: c+2.2 
0.0_2 


4.4«;« 

o 7+1.0 

o 1+0.9 
■'•-'--0.8 

5-oi};? 


l-849iH'oI(6.71) 
1.848l°:»^(6.71) 
1.858lH?2(6.67) 
l-857tHoi(6-68) 
1.847iHo2(6.71) 


< 16.4(< 2659) 

< 10.0(< 1623) 

7«(117ltf7f) 
< 4.6(< 746) 


n 1+2.3 

■3-l_1.4 

^•"-2.0 
4.5t^i 


3.0+}-' 

2.8+}- 

4.1+}-' 

1 9+°-' 
i.y_Q , 


Note. — 


EW is the equivalent 


width, and the line flux is in 


units of 10 


ph cm~^ . 


Errors are quoted at 90% confidence. 



Table 4. Best-fit parameters from xstar photoionization models of the near-dip and 

far-dip phases. 



phase 


log{C) 
(ergs cm s" 


log(AfH) 
-1) (cm-2) 


V 

(km/s) 


log(0 log(ArH) 
(ergs cm s"-"-) (cm~^) 


V 

(km/s) 


near-dip 
far-dip 


4 9+0.8 


Hot Absorption 
23 4+°-'5 
90 Q+0.3 
^■-'•■-'-0.3 


ooc: + 346 
■-''-'■'-318 
9CO+306 
^"■'-310 


Warm Absorption 
0-1+0.2 990-1-0.3 
■'•-'-0.2 ^^"'-0.4 
4+0.2 99 9-1-0.2 
■'•^-0.1 ^^-^-0.2 


-2761^71 
1028t^890 



Note. — ^ is the photoionization parameter from XSTAR. 'Wh" is the equivalent Hydrogen column density in the associated 
ionized region, "w" is the measured velocity shift, (-) blueshift towards us and (-|-) redshift away us. Errors are quoted at 90% 
confidence. 
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9.11]'^Q24 X 10^^ cm~^ and T = 2.261;qo5. Figure [6] shows that the equivalent hydrogen column 
density of the local cold absorber varies significantly over the orbital phase, while the partial 
covering factor remains relatively steady. This variation of the hydrogen column density 
of the local cold absorber is very similar to the phase-dependent variations seen for the 
absorbed power-law fiux in the alternative blackbody-|-power-law models described above; 
therefore, we also use a sine function to fit the points of Fig. [6]-(III). At 90% confidence, we 
derive a period T = 47 ± 7 ks, similar to the fits above. Again, we are seeing evidence that 
the observed orbital modulations are predominantly driven by absorption variations, rather 
than by intrinsic source variations, regardless of whether we fit a blackbody plus power-law 
model, or instead fit a partial covering model of a single component (i.e., the power-law). 

We also investigate Fe XXV and Fe XXVI absorption line changes over these smaller 
time segments by fitting Gaussian lines to the absorption features while fixing the continuum 
parameters to the phase-specific values of Table [51 (Here, the HEG spectrum is binned to 
0.005 A bin-^) Table M and Figs. H and [5]-(IV) show that the Fe XXVI velocity width is 
narrowest during phases c and d that occur immediately after dipping events. In contrast, it 
appears that it is the Fe XXV fiux and velocity shift that change most dramatically over the 
orbital period, supporting our interpretation that these lines originate in different regions 
(see f|3]). This result is consistent with the findings from previous sections (see, e.g.. Table H]) 
whereby a comparison between xstar modelling of the the near-dip and far-dip regions show 
that the largest changes come from the less ionized "warm" region of log^ ~ 3. 



3. Discussion 

Our Chandra-HETGS spectral study of 4U 1624—490 over its orbital period confirms 
the detection of H- and He-like Fe XXV and Fe XXVI absorption. Detailed photoionization 
modelling with xstar predicts additional lines from Ca XX, Ar XVII, and S XVI. Of partic- 
ular note is not only the confirmation of these lines with detailed photoionization modelling, 
but also that the variability analysis strongly suggests that these lines can be identified as 
originating from the accretion disk corona (dominated by high ionization plasma giving rise 
predominantly to Fe XXVI) and the outer accretion disk (as indicated by the fiux evolution 
of the lower ionization variable Fe XXV absorption). 

While the time-averaged spectrum over the 4U 1624—490 orbital period can be modeled 
reasonably well with one ionization zone, variability studies instead suggest that a model with 
two ionization zones is a better description of the data. This interpretation has important 
implications for geometric scenarios and our ability to m ap the emitt i ng reg ions as a function 



of orbital phase. In an earlier analysis of this data set, llaria et al.l (120071 ) claimed velocity 
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Table 5. Best-fit parameters over ~10 ks time step. 





kT 


uFlux 


Flux 


r 


uFlux 


Flux 


E 


Width 


Flux 


phase (10^^ cm^^) 


(keV) 












(KeV) 

(A) 


(KeV) 
(10^ km/s) 





blackbody 



power-law 



Fe emission feature 



a 


7. 


fj+0.7 


1, 


,35«: 


,09 
,05 


11, 


.7tl' 


1 

.2 


6, 




,5 
.3 


0, 


O + O, 

■S-o 


7 
.8 


7 7+4.5 
' ■ '-3.4 


4, 


.9tS: 


,3 
,8 


a c-T+0.12 

(1.858) 


0.25lg:?? 

(His) 


1, 


7+0. 
■ '-0, 


b 


7, 


o+l.O 
•^-0.3 


1, 


,32±°: 


,07 
,06 


11, 


.3tl' 


,6 
.0 


6, 


n + O, 


,9 
.8 


0, 




,9 
.5 




4, 


.5t?: 


,0 
,1 


"■^^-0.25 

(1.925) 


< 0.8 

« 37) 


0, 




c 


7, 


4+0.4 
■^-0.3 


1, 


,46«: 


12 
,11 


12, 


.Oil 


,4 
.2 


6, 


,811: 


,3 
.4 


1, 




,3 
.0 




2, 




,6 
,6 


(1.781) 


n iq+o.6o 

(49ilf ) 


0, 




d 


8, 


Q+1.4 


1, 


,481°: 


12 

,11 


12, 




,0 
.2 


6, 




,5 
.3 


1, 




,1 
.4 


« (;+l3.9 

»-0-5.g 


3, 


5+1, 

■0-0, 


5 

,8 


(1.794) 


o.3otO:i| 

(13t«6) 


1, 


■0-0, 


e 


7, 


4+0.6 


1, 


,321°: 


,08 
,06 


11, 


9+2, 


,9 
.8 


6, 


1+0, 

■-'--1 


,7 
.5 


0, 


7+0, 
■ ' -0 


,7 
.7 




5, 


2+1, 


,6 
,8 


6.58«:1« 
(1.884) 


0.211°:?? 


1, 


Q + 0. 


f 


7, 




1, 


,3i«: 


14 
,05 


11, 




,3 
.6 


6, 


O + O, 
■•^-0 


,7 
.9 


1, 




,1 
.5 




3, 


o + l, 

■«-a. 


1 

,8 


D.Dy„Q 22 

(1.853) 


0.28i°:?^ 
(13tf) 


1, 


■otS: 


g 


7, 


7+0.7 
'-0.5 


1, 


,26±«: 


,06 
,05 


11, 


.8tl 


,7 
.6 


6, 


9+0, 
■^-1 


,9 
.0 


1, 


.It?' 


,8 
.1 


5 ^+*-5 
■J-'J-3.1 


3, 


2+1, 


,0 
,1 


6.82«:?1 
(1.818) 


0.26t°:?^ 


1, 


■4iJ: 



Note. — A^H is the LOS cold Hydrogen column density. uFlux and Flux refer respectively to the unabsorbed and absorbed 
integrated flux [10"^" ergs cm"-^ s"^] over 0.6-10.0 keV. lunc is the flux [10^^ ph cm^^ s^^] of the broad Fe emission line. 
Errors are quoted at 90% confidence. 



Table 6. Best-Fit parameters of Fe H-Iike and He-like absorption lines. 





Wavelength 


Width 


EW 


Flux 


Wavelength 


Width 


EW 


phase 


(A)(keV) 


(mA)(km/s) 


(mA) 




(A)(keV) 


(mA)(km/s) 


(mA) 



Flux 



Fe XXVI 



Fe XXV 



1.780lo 
1.780lg: 
1.782lg 
1.775lg: 
1.776111! 
1.7771J 
1.780lo 



■r2(6.97) 
.r3(6^97) 
•°o^(6.96) 
°°3(6.99) 
■gg5(6.98) 
°J?(6.98) 



< 7.4(< 1250) 

< 6.2(< 1042) 

< 4.4(< 749) 

< 3.9(< 667) 
8tl{l377tlll) 
< 15.0(< 2532) 

9.6t^;7(2028t?i?i) 



4. it? 



7 



,+ 1.9 

'-2.5 

14^5il| 



n 

4.2t? 
3.9i?. 

6.7l2' 
7.11^ 



3.2^ 



1.845tHS2(6.72) 
1^852iHJ6(6^69) 
1.850tE;:0?o(6-70) 
1.860iHo4(6.67) 
l-845lE!:°S|(6.72) 
l-845i°:°S5(6.72) 



8.3tl:l 1.850iH?,^5(6.70) 



< 6.2(< 1006) 

< 20.6(< 3342) 

13l^7{2077l?J«? 

< 6.9(< 1110) 

< 10.0(< 1626) 

< 4.9(< 803) 

< 7.5(< 1218) 



4.7 
6.0 
13.1 
5.7 
< 6.5 
3.3 
3.8 



+ 1.8 
2.1 
4.4 
3.4 
+ 7.0 
5.6 
+ 1.9 
2.2 



2.2 
2.4 
+2.2 
2.7 



3.2 
3.8 
7.1 
3.8 
< 4.2 
2.0 
2.2 



+1.2 
1.4 
+2.8 
2.2 
+3.8 
3.0 
+ 1.2 
1.4 



+ 1.3 
1.5 
+ 1.3 
1.5 



Note. — EW is the equivalent width and the unit of flux is 10 * ph cm ^ s ^ . Errors are quoted at 90% confldence. 
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1.6 



Phase (g) 

1.7 1.8 ^ 1.9 

Wavelength (A) 



Fig. 4. — HEG spectra (black) and best-fit model (red) showing the evolution of the Fe XXV 
and Fe XXVI absorption lines over the 4U 1624—490 orbital period. As in Fig. [3], the vertical 
green lines indicate the laboratory wavelengths of Fe xxv (1.850A) and Fe xxvi (1.778A). 



- 19 - 



3 
a 

G 
3 

O 
U 



00 



X 



a 

o 
si 



Light Curve of 4U 1624-490 

Ma 



t\J\ 



o 1 



o 

00 



in 
I 



o 
in 



in 
o 



(I) 

f I g 



unabsorbed blackbody 
unabsorbed poweiiaw 
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40 
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Fig. 5. — (I) Tracking the 4U 1624—490 hght curve over 10 ks time steps, shows that 
the unabsorbed fluxes of both the blackbody and power- law components (II), and to a 
lesser extent the absorbed flux of the blackbody component (III), remain relatively steady, 
whereas the absorbed power-law component (III) reaches a minimum during the periods 
near dipping events. Both components can be described by sine waves (blue curves). (IV) 
Intensity evolution of the Fe absorption lines. The blue curve is also a sine wave fltting 
to the intensity of Fe XXV absorption line. (V) The evolution of the intensity ratio of the 
absorption lines shows a quasi-sinusoidal variation. In this flgure only, all error bars shown 
are at the la confldence level. The vertical bars in subflgure (I) indicate the different phases 
of flux evolution. 
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Light Curve of 4U 1 624-490 (I) 




20 40 60 

TimeMJD(10** s) 



Fig. 6. — The evolution of (II) the partial covering factor and (III) the equivalent hydrogen 
column density of the local absorber, which reaches a maximum during the periods near 
dipping events. Similar to the fits in Fig[5|, the equivalent hydrogen column density also can 
be described by a sine wave with a period of T = 47 ± 7 ks (dashed curve). In this figure, 
all error bars shown are at the Icr confidence level. 
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widths of up to 3500 km s~^, to which they attribute an origin between the external region 
of the ADC and the disk edge. In this interpretation, the broadening would be due to bulk 
motion or turbulence connected to the coronal activity. Our variability results indicate that 
neither bulk motion nor turbulence is necessary to produce the line widths. We arrive at this 
conclusion based on variability studies where line velocity widths are seen to change over 
relatively short time scales, and that in some phases the velocity widths are comparable to 
the HEG FWHM resolution (i.e., Fe xxvi during the near-dip phase, and Fe xxv during the 
far-dip phase) and therefore are easily described by thermal broadening. 

To take our two-region model further, we consider it in the context of simple arguments 
relating the ionization parameter to the location of the absorbing plasma. As pointed out by 
Tarter. Tucker fc Salpeterl (119691 ). the ionization parameter ^ can be defined by the relation 
^ = Lxn~^R~^, where Lx is the total un-absorbed source luminosity (= 4.9 x 10^^ ergs s~^ 
between 1-10 keV during the epoch of our 4U 1624—490 observation) and R is the distance 
to the absorbing plasma. The electron particle density Uc = Ni{/AR is characterized by the 
thickness of the absorbing medium, AR, and the equivalent Hydrogen column of the ionized 
region. Using this, we consider the possible range of R for the plasma giving rise to the 
observed ionized absorption lines. 

A sketch of a scenario explaining the major results of such an analysis is shown in 
Fig. [71 As discussed above, the observed orbital variability of the different spectral com- 
ponents points towards an exte nde d ADC similar to the emission geo metry proposed by 
White. Stellar k Parmail Jl988h or Ichurch fc Balucihska- Church! J^OoJ (see Ichurchi 12004 
for a summary). In this model, the ADC is large and therefore only slight variations of 
its components with orbital phase are expected, in line with our observation that the hot 
plasma component is constant over the orbit. This location of the hot component is con- 
sistent with the spectral analysis: Our modelling efforts show that the "hot" ga.s compo- 
nent can be described by ^hot ~ 10^'^ ergs cm s~^ (T ^ 3 x 10^ K) ionizing a plasma of 
column ^ 2 X lO^^cm"^. "Naively" assuming a geometrically thick AR ^ R, a distance 
i? ~ 3 X 10^° cm is obtained, which is comparable to the size of the Accretion Disk Corori a 
(ADC) derived from the relation i?ADC = L°^^^^°-^^ of IChurch &: Balucihska- ChurchI J2004h . 
A similar value for an extended ADC was recently determined by Schulz et al. (j2008b ) while 
measuring Doppler broadening in several broad emission lines in Cygnus X-2. 



In contrast, ^warm changes more dramatically over the 4U 1624—490 orbital phase with 
the largest differences seen in a comparison between the far-dip (^ ~ lO^ '^ergscms"^; T 
1.0 X 10^ K) and near-dip (^ ~ 10^'^ ergs cm s~^; T ^ 0.9 x 10^ K) phases. This component 
cannot originate in the ADC: If one assumes AR ~ R and a column density of iVn ~ 2 x 
lO^^cm"^ as appropriate for the "warm "component (Table H]), the data require R ~ 10^^ cm. 
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which is significantly larger than the likely accretion disk radius of ~ 10^^ cm. If, however, a 
geometrically thin plasma of AR/R ~ 0.1 is assumed for this "ti;arm" component, absorber 
distances between 1.8 x 10^^ cm (near-dip) and ~ 1.1 x 10^^ cm (far-dip) are obtained. Since 
the radius of the accretion disk is ~ 1.1 x 10^^ cm as obtained from the Roche Lobe size of 



the neutron star (jChurch &: Balucihska-Churchl |2004| ). it is therefore likely that the warm 



component originates in the accretion disk rim, e.g., where the accretion stream impacts 
onto the disk and where gas is irradiated by the neutron star. 

We note that this interpretation of the absorbers is also consistent with the varia- 
tion of the broad-band s pectra l par ameters shown in Fig. [5l In terms of the model of 



White. Stellar fc Parmarl (119881 ) and IChurch fc Balucihska-Churchl (120041 ). the blackbody 



component is taken to originate from the neutron star, while the power-law component is 
interpreted as being due to Comptonization in the extended ADC. In this interpretation the 
power-law component is expected to exhibit a quasi-sinusoidal variation and to be weakest 
during the times of dipping when the ADC is partly covered by the warm absorbing material. 
A consequence of this model is also that the ratio Fe xxv/Fe xxvi shows a maximum during 
the dipping, since the optical depth in Fe xxv is largest during that time. 

We also note that the evolution of the spectrum can be modeled with partial covering 
of a single power-law component, with the local absorber varying over orbital phase. In such 
a model, the soft-end of the uncovered power-law replaces the blackbody component. This 
offers a slightly different explanation of the relative strength of the variation of the absorbed 
power-law component, relative to the blackbody. The "covered power-law" is that portion 
of the ADC at low latitudes, where obscuration by the disk is expected to be greatest, and 
the ADC is expected to be widest. The "uncovered power- law" arises from a geometrically 
narrower region at high latitudes, where obscuration by the disk is the least. For either the 
blackbody plus power-law scenario, or the partial covering scenario, our analysis shows that 
variations are predominantly driven by changes in obscuration, rather than any intrinsic 
variation of the components. 



4. Summary 

We summarize the major conclusions and discoveries of this paper: 

• A re-analysis of the Chandra -HFiTGS spectrum of 4U 1624—490 based on photoion- 
ization modelling and variability studies, shows that while the continuum spectrum 
can be modeled using a blackbody plus power-law, a slightly more preferred model is 
a single F = 2.25 power-law partially (71%) covered by a local absorber of column 
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Fig. 7. — The 4U 1624—490 geometry in the context of our hne-of-sight view as determined 
from absorption hne variabihty studies. Our studies indicate that the hotter component 
dominates the ADC, whereas the cooler gas is associated with the outer disk rim, which 
intersects our hne-of-sight. a-g refer to the 4U 1624—490 orbital phase as indicated by the 
light curve of Fig. [1^. 



-24- 



density A^h, Local = 8.11q;q x lO^^cm"^. The total 1.6 x 10^^ cm"^ large line-of-siglit 
(^H, Local + ^H, Gal) columii therefore can attenuate the power-law up to 5 keV, even 
though the X-ray spectrum can be seen down to 1.7 keV. For spectra with limited 
(e.g. 0.5-10 keV) band coverage, the severe absorption at the low energies coupled 
with limited counts at the high energies can flatten a F = 2.25 power-law to ap- 
pear as F ~ 1.2. The source luminosity at the epoch of our observations is esti- 
mated to be 4.9 x 1 ^^ erg s s~^ for a distance of 15 kpc as recently determined by 



Xiang. Lee fc NowakI (120071 ). 



• We report the discovery of a possible quasi-sinusoidal modulation with period P = 
43l9^ o^^^ the ~ 76 ks orbital period, with the caveat that we have observed 
4U 1624—490 only over one binary orbit. There is strong indication that variations 
in observed spectral properties are predominantly driven by changes in obscuration, 
rather than any intrinsic variation of the power-law or black-body components. 

• We confirm the detection of i onized absorption from Fe XXVI, Fe XXV, and Ca XX, 



which have been detected by llaria et al.l (120071 ). Detailed photoionization modeling 



with XSTAR has allowed us to derive the strengths of other absorption lines of Mg, Si, 
S and Ar. 

• Our detailed variability studies of detected strong absorption lines point to a two- 
temperature plasma for their origin: a highly ionized component of ^hot ~ 10^ '^ ergs cm s~^ 
(T~ 3.0 X 10^ K) associated with an extended accretion disk corona of i? ~ 3 x 10^*^ cm, 
and a less ionized more variable component of ^ ~ 10^'^ ergs cm s~^ (T~ 1.0 x 10^ K; 
far-dip) and ^ ~ 10^'^ ergs cm s^^ (T~ 9.0 x 10^ K; near-dip) coincident with the ac- 
cretion disk rim. This model is offered in favor of bulk Comptonization, as previously 
suggested, as an explanation for observed line broadening. 

• We also report Fe XXV and Fe XXVI absorption line variations (in wavelength, strength, 
width) over the 4U 1624—490 orbital period. This, with previously noted variations 
have allowed us to map plasma properties to our viewing of 4U 1624—490 over its 
~ 76 ks orbital period (see Fig. [7]). 
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